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Abstract 

We report the first observation of a number of decay modes of the B meson, namely B° — > D*~ (5tt) + , 
B + — > D*~(4tt) ++ and B + — > D (5tt) + , where (nir) implies the combination of n charged pions. The 
analysis is based on a 140 fb _1 data sample collected at the T(4S I ) resonance with the Belle detector at 
KEKB. We measure B{B° -► D*- (5tt)+) = (4.72±0.59±0.71) x 1(T 3 , B(B + -► D*-{Air) ++ ) = (2.56±0.26± 
0.33) x 10~ 3 and B(B + — ► D*° (57r) + ) = (5.67±0.91±0.85) x 10~ 3 . We also provide improved measurements 
of the branching fractions for the decay modes B° -► D*~ (3vr) + , B + -»• D*° (3vr)+ and B° -> D*° (4vr)°. 

PACS numbers: 13.25.Hw 
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A large fraction (~ 35%) of B meson decays is due to decay modes which are still unknown. While 
a sizable fraction of B's decay into semileptonic final states (~ 24%) decays into hadronic final 
states are dominant. At present, however, only half of these final states correspond to measured 
exclusive hadronic decays. 

On average, B mesons decay into a large number of particles: the mean charged particle multi- 
plicity in hadronic B decay is measured to be 5.8 ± 0.1 j2j. Decay products of D mesons make a 
significant contribution, but B decays to charmed states with a large number of accompanying pions, 

U\ 1—1 

B — ► D (nir), where the 7r are charged pions and n = 3,4, are also known [3j|. The invariant mass 
distribution of the multi-pion system, and the resonant decomposition for such decays, are important 
for the study of factorization |4|. 

In this paper, we present a study of inclusive B — ► D (rwr) final states, where n = 3,4 and 5, 
and measure branching fractions for six decay modes. Inclusion of charge conjugate modes is implied 
throughout this paper. 

The analysis is based on a 140 fb" 1 data sample at the T(AS) resonance (10.58 GeV) and a 16 
fb" 1 data sample 60 MeV below the T(4S) peak (referred to as off-resonance data), collected by the 
Belle detector |H[ at the asymmetric e + e~ collider KEKB The data sample contains 152 million 
BB events. 

The Belle detector is a general purpose magnetic spectrometer with a 1.5 Tesla magnetic field 
provided by a superconducting solenoid. Charged particles are measured using a 50 layer Central 
Drift Chamber (CDC) and a 3 layer double sided Silicon Vertex Detector (SVD). Photons are de- 
tected in an electromagnetic calorimeter (ECL) consisting of 8736 CsI(Tl) crystals. Exploiting the 
information acquired from an array of 128 time-of-flight counters (TOF), an array of 1188 silica 
aerogel Cerenkov threshold counters (ACC) and c?£'/(ix-measurements in CDC we derive particle 
identification likelihoods £ n /K- A kaon candidate is identified by a requirement on the likelihood 
ratio Ck/(C-k + Ar) such that the average kaon identification efficiency is ~ 90% and the pion fake 
rate is ~ 9%. Similarly charged pions are selected with an efficiency of ~ 91% and the kaon fake 
rate is ~ 10%. We select charged pions and kaons that originate from the region |Ar| < 0.2 cm and 
\Az\ < 4 cm with respect to the run dependent interaction point, where Ar, Az are the distances 
of closest approach of tt/K tracks to the interaction vertex in the plane perpendicular to the beam 
axis and along the beam axis, respectively. All tracks compatible with the electron hypothesis (~ 
0.2% fake rates from pion/kaon) are eliminated. No attempt has been made to identify muons, which 
represent a background of about 2.7% to the pion tracks. Candidate 7r° mesons are identified as a 
pair of isolated ECL clusters with invariant mass in the window 118 MeV/c 2 < M 77 < 150 MeV/c 2 . 
The energy of each photon is required to be greater than 30 MeV in the barrel region, defined as 
32° < 6* 7 < 128°, and greater than 50 MeV in the endcap regions, defined as 17° < 6* 7 < 32° or 
128° < 9 1 < 150°, where 9 1 denotes the polar angle of the photon with respect to the direction 
opposite to the e + beam. A mass constrained fit is applied to obtain the 4-momenta of 7r°'s. 

Beam gas events are rejected using the requirements \P Z \ < 2 GeV/c and 0.5 < E w - ls j*Js < 1.25, in 
the T(45') rest frame, where P z and E vis are the sum of the longitudinal momentum and the energy 
of all reconstructed particles, respectively, and is the sum of the beam energies in the T(4S") rest 
frame. 

The D meson is reconstructed through its decay to K + tt~ . A vertex constrained fit is performed 
and the invariant mass is required to be within 17 MeV/c 2 (~ 3.5 cr) of the nominal D mass. D*~ J s 
are then reconstructed by combining the D° with a slow charged pion with | Ar| < 2.0 cm and | Az\ < 
10.0cm with respect to the D vertex (D vertex resolutions are a T ~ 0.026 cm and a z ~ 0.016 cm). 
The signals of B — > D (nn) are reconstructed through the decay chain D —> D tt . Candidate 
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D*~ (D* ) are selected when the reconstructed mass difference between D* and D is within 2.0 (2.4) 
MeV/c 2 of the nominal mass difference AM, which corresponds to ~ 3.5 (3.0) a resolution on AM. 
A kinematic fit with the nominal D mass is applied to obtain the 4-momenta of the D candidate. 

The D candidate is then combined with n pions to reconstruct the B meson. All n pions are 
required to satisfy |Ar| < 0.2 cm and |Az| < 0.8 cm with respect to the D vertex. Continuum 
(e + e~ — > qq, where q = u,d,s,c) events are suppressed with the criterion, | cos thrust I < 0.8, where 
^thrust is the angle between the thrust axis of the B candidate daughters and the thrust axis of the 
remaining tracks and isolated ECL clusters. 

The signal can then be identified by two kinematic variables calculated in the T(4S) rest frame. 
The first is the energy difference, AE = Ed* + 127=i — -^beam, where Ed* is the energy of 
the D candidate, E ni is the energy of the zth pion in the rm system and -E be am = \fs/^ (Fig. 

[TJ). The second variable is the beam-energy constrained mass, M bc = J E^ cam — \Pd* + 127=1 l 2 > 
where Pd* and P 7Ti are momentum vectors of the D* candidate, and the ith pion in the nir system. 
Typical resolutions for these variables are 8.5 MeV and 2.7MeV/c 2 , respectively. In the extraction 
of the signal yield, we require 5.273 GeV/c 2 < M bc < 5.288 GeV/c 2 and M bc > 5.27 GeV/c 2 
for D*~ (nn) and D (nn), respectively, and we fit the AE distribution from —150 MeV to 150 
MeV. In many B decay analyses, the AE distribution includes peaks or other structures due to 
related B decays with an additional particle, one particle less than in the mode under study or 
misidentification of a particle from a topologically similar decay mode. In this case, we do not 
observe such structures within our fitting range. We do not use the M bc distribution to obtain 
signal yields, because peaking backgrounds in that distribution cannot be distinguished from signal. 
Selected events contain multiple B candidates with a multiplicity depending on the signal channels, 
which varies from 1.2 to 1.8. For each event we choose a unique B candidate, taking the combination 
resulting in the minimum value of ((M bc - M B )/(TM hc ) 2 + (AM d /ct Md ) 2 + (AMd* /(Tm d *) 2 - 

We have studied continuum events and other B decays as possible sources of background. Back- 
ground due to continuum events is studied by analysing the 16 fb _1 off-resonance data and there we 
do not find any peak near AE — 0. 

B decays (both signal and background) are studied with Monte Carlo (MC) event samples. MC 
events are generated using the QQ event generator 0] with a phase space distribution for the tvk 
system, and the response of the Belle detector is simulated by a GEANT3-based program [§| . The 
simulated events are then reconstructed and analysed with the same procedure as is used for the real 
data. 

We have investigated the possibility of reconstructing B — > D* (n'7r)(mir°) channels as 
D*~(°)(n7r) due to the loss and/or addition of pions, where n = 3, 4 or 5, n' = n, n ± 1 and m = 
or 1 and here we observe a linear background without any structure around AE = 0. We have also 
studied B — > D* (— > D 7)(n7r) MC events reconstructed as B — > D* (— ► D ir )(nir) and D*~^ (nir) 
events reconstructed as D* ^ (nir). In each case, we find a AE distribution peaked at zero, but 
with a larger width than the corresponding signal distribution. We also consider backgrounds due 
to B — > D {rm) and B — ► _D°(mr)7r° with n > 4. Such decays have not yet been observed; assuming 

branching fractions equal to those of the corresponding B —> D (nir) and B — > D ((n+l)7r) modes, 
respectively, we find a small background to D*~(mv) for n = 3,4, and a negligible contribution to 
the other final states. Finally, we fit MC signal distributions without and with contributions from 
the various feed-across backgrounds. The ratio of those two signal yields (F r ) depends on the signal 
channels and varies from 0.99 to 0.87. Observed signal yields are corrected with the corresponding 
F r when extracting branching fractions. 
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D*° (3tt) 



We have studied the following six B decay modes: B° — > D* (37r) + , 
B+ _> D*-(47r)++, 5° -> D*° (47r)°, 5° -»■ D*~ (5tt)+ and 5+ -> 77*° (5tt) + . Figure □ shows the AE 
distributions for these decay modes. Statistically significant structures near AE =0 are observed, 
We have also checked the corresponding AE distributions for events in the M^ c sideband region 
(5.23 GeV/c 2 < M bc < 5.26 GeV/c 2 ): no structure is observed. For all the decay modes under 
study, backgrounds are fitted with a linear function. The signal shape is modelled with a sum of two 
Gaussian distributions for B — > D* (rwr). The D (rwr) signal is fitted with a sum of the Crystal 
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and a Gaussian. The n and a parameters of the CB and the fractional area of the second Gaussian 
are fixed from fits to the MC sample. Signal yields obtained from the fit are summarised in Table |U 



TABLE I: Measured signal yields and branching fractions. 



Channel 




$nal yield 




B x 10 3 


B°- 


-> D*~ (3vr)+ 


2554.9 


± 


58.4 


6.81 


± 


0.23 


± 


0.72 


B+ 


-> D* (4ir)++ 


445.9 


± 


31.7 


2.56 


± 


0.26 


± 


0.33 


B°- 


-> D*~ (5vr)+ 


449.6 


± 


32.0 


4.72 


± 


0.59 


± 


0.71 


B+ 


-> D*° (3vr)+ 


1728.4 


± 


52.9 


10.55 


± 


0.47 


± 


1.29 


B°- 


-> (4vr)° 


280.6 


± 


36.5 


2.60 


± 


0.47 


± 


0.37 


B+ 


- D*° (5tt)+ 


341.2 


± 


40.6 


5.67 


± 


0.91 


± 


0.85 



The signal efficiency depends on the invariant mass (M nn ) of the me system. Signal efficiencies 
are calculated from MC event samples in 100 MeV/c 2 bins of M n7r . 

In order to obtain the signal yield in each M n7T bin, we perform a fit to AE distributions in each bin 
to avoid the possible uncertainty due to different background shapes in signal and sideband regions 
of Mb c . Efficiency corrected M nn spectra are shown in Fig. |21 Invariant mass distributions of the 37r 
system (M( 37r )) in B —>■ D (3%) decays show clear evidence for the presence of ai (Fig. Efa,d)). The 
M 47r distributions in Fig. |2(b,e) do not show any resonant structure. M 57r distributions in Fig. |2fc,f) 
show a peak around 2 GeV/c 2 which, however, does not correspond to any known resonance. We 
also search for narrow resonances in the M 5n distributions, after the subtraction of AE sidebands. 
Apart from a peak due to the D+, no narrow structure is observed. 

Inclusive branching fractions obtained for the six modes are listed in TableHJ where the first error is 
statistical and the second error is systematic. This is the first measurement of the branching fractions 
of B + — > D*~(47r) ++ , B° — > D*~ (57r) + and B + — > T>*° (5n) + . The branching fractions for the decay 
modes B° -> D*~ (3tt) + , B+ -> D (3vr)+ and B° -> D (4tt)° are measured with better precision 
than in previous studies |lj. We do not exclude contributions from exclusive channels B — > D*D+, 
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FIG. 1: AE distributions for six D*(nvr) combinations: (a) B° -> D*~ (3vr) + , (b) 5+ -> L>*-(47r)++, (c) 
5° -> L>*- (5vr)+, (d) -> (3vr)+, (e) B° -> (4vr)° and (f) B+ -» (5vr)+. Points with error 
bars are the observed events in data, solid lines are the results from the fit and dashed lines represent the 
background components. 



B — > D* D* + , B — > D* D + , B — > D*D* + Ks etc. The largest contributions of this type to inclusive 
D*(3tt) and D*(5vr) are expected to come from D*D+, with D+ -> (3vr)+ (B = 1.01 ± 0.28% Q) 
and £>+ — > (5tt) + (S = 0.65 ± 0.18% [l|), respectively. We conduct a search for these modes by 
fitting invariant mass distributions of nir, M n7T for signal events (|AJ5| < 30 MeV) and find yields of 
27.8±7.7, 15.6±5.4, 8.4±3.9 and 11.3±5.4 in B° -> D*~ (3tt) + , 5+ -> D*° (37r) + , 5° -> £>— (5tt) + 

and _B + — > D (5tt) + channels, respectively. Within the statistical errors, these results are consistent 
with PDG expectations. The expected contribution from other modes are smaller than B — > D 
by a factor ~ 4. 

The systematic uncertainty is obtained from a quadratic sum of nine terms, which are shown in 
Table |H] the uncertainty in (a) the track finding efficiency, ranging from 1% for high momentum 

tracks to 8% for 80 MeV/c pions, calculated from partially reconstructed D*~ — > D°(^ Ks(~^ 
7r + 7r~)7r + 7r~)7r~ events and a track embedding study; (b) the slow tt finding efficiency; (c) K/tt 

selection efficiencies (PID), calculated using events; (d) branching fractions 

of D*~(D* ) -> D°7r-(°), D -> TT+tt- and tt — > 77; (e) the uncertainty in the feed-across from other 
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FIG. 2: Efficiency corrected nir invariant mass spectra for (a) B° — > D* (37r) + , (b) B + — > D* (47r) ++ , (c) 
5° -► L>*- (5vr)+, (d) 5+ (3tt)+, (e) B° -► (4tt)° and (f) 5+ -► (5tt)+. 



TABLE II: Contributions of systematic uncertainties (in %). 





D*- 


ZP° 




4vr 


57T 


3vr 


4vr 


5vr 


Track finding 


8.7 


11.0 


12.7 


5.4 


6.6 


7.9 


Slow 7T° 








7.0 


7.0 


7.0 


PID 


4.9 


5.0 


4.8 


5.0 


5.2 


5.3 


Branchings 


2.5 


2.5 


2.5 


5.3 


5.3 


5.3 


Feed- across 


0.8 


1.0 


2.0 


0.9 


1.4 


2.2 


Fitting 


0.5 


1.2 


1.4 


2.4 


4.0 


5.5 


Selection 


2.5 


3.7 


4.9 


3.3 


5.7 


4.7 


MC 


0.8 


0.8 


2.2 


1.2 


1.1 


1.7 




0.5 


0.5 


0.5 


0.5 


0.5 


0.5 


Total 


10.6 


13.0 


15.0 


12.3 


14.1 


15.0 
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decay modes which is calculated by changing the relative branching fractions of backgrounds to signal 
by one sigma of their errors; (f) the uncertainties due to the choice of signal shape and histogram 
binning. AE distributions are also fitted with (i) single Gaussian (CB), (ii) modified Gaussian (CB), 
where Gaussian power has changed from 2 to 1.4, (Hi) asymmetric Gaussian, (iv) modified Gaussian 
(CB) + Gaussian and (v) asymmetric Gaussian + Gaussian for D*~ (rar) (D (nir)) signals. Ratio 
of fitted yields in data and MC signals are calculated for these functions and the variation of these 
ratio is taken as the systematic error due to fitting function; (g) the uncertainties due to selection 
criteria, which are calculated by varying those criteria by one sigma of their errors, (h) limited MC 
statistics and (i) the uncertainty in the total number of BB events (N BB -). 

In summary, we have made first observations of decay channels B° — ► D*~ (57r) + , 
B + _> D*-(4n) ++ and B + -> D (5tt)+ usine; 152 million BB events. We measure inclusive branch- 
ing fractions for these three decay modes. We have also made precise measurements of the branching 
fractions for the decay channels B° -> D*~ (3tt)+, B+ -> D*° (3tt)+ and B° -»■ D*° (4vr)°. 
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